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Nucleon electromagnetic form factors from lattice QCD 
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Abstract. It is imperative that lattice QCD serve to develop our understanding of hadron structure and, 
where possible, to guide the interpretation of experimental data. There is now a great deal of effort directed 
at the calculation of the electroweak form factors of the nucleon, where for example, measurements at 
Jefferson Laboratory have recently revealed surprising behaviour in the ratio Ge/Gm. While, for the 
present, calculations within the framework of lattice QCD are limited to relatively large quark mass, there 
has been considerable progress in our understanding of how to extrapolate to the chiral limit. Here we 
report the results of the application of these techniques to the most recent form factor data from the 
QCDSF Collaboration. The level of agreement with all of the form factors, for below 1 GeV^, is already 
impressive. 



PACS. 12.38.Gc Lattice QCD calculations - 11.30.Rd Chiral symmetries 
form factors - 14.20.Dh Protons and neutrons 
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1 Introduction 

The electromagnetic form factors of the nucleon provide a 
fundamental constraint for any theoretical description of 
the structure of the nucleon p. Even though they have 
been studied experimentally for more than 50 years, recent 
experiments at Jefferson Lab j2] have revealed surprising 
new behaviour in the ratio Ge/Gm for the proton. In 
addition, there are a range of other new results on the 
neutron electric and magnetic form factors IT - this is an 
exciting and rapidly developing field. 

In this context it is vital that lattice QCD, our only rig- 
orous method of solving non-perturbative QCD, be used 
to inform our understanding of this data and the vari- 
ous models of hadron structure that are used to describe 
it. There was a relatively long hiatus in calculations of 
nucleon form factors in lattice QCD after the pioneering 
work of Leinweber and collaborators in the early 90's ^ 
I3n|- However, the activity has intensified in the last few 
years [3|S1|H] . In spite of this activity, limitations in com- 
puter speed mean that we are currently limited to lattice 
simulations at quark masses a factor 5-10 higher than the 
physical values and these calculations are currently made 
in quenched approximation (QQCD). 

If one is to compare these state of the art simulations 
with experiment, it is necessary to make an extrapolation 
as a function of quark mass to the physical region [lUL 
lll| . This extrapolation is made non-trivial by the non- 
analytic behaviour as a function of niq which follows from 
the fact that chiral symmetry is dynamically broken in 
QCD. Until recently the absence of data has meant that 



efforts at chiral extrapolation of form factors has been 
focussed on baryon magnetic moments jl2lll3lll4lll5| and 
charge radii p.6 ,17 . The focus on non-analyticity jJHI of 
hadron electromagnetic properties has inspired consider- 
able investigation of the formal constraints in both full 

The practical issue is then how to incorporate these 
formal chiral constraints into a practical chiral extrapo- 
lation, given that the radius of convergence of the for- 
mal expansion dictated by chiral perturbation theory is 
rather small 24 . In the case of the nucleon mass, where 
there is extremely accurate data in full QCD from the CP- 
PACS Collaboration [221, one can formally demonstrate 
the model independence of the extrapolation procedure 
•24"26'"27' . For the magnetic moments the data is only now 
improving to the point where the model independence of 
the choice of finite range regulator |2H1 can be examined 
|12| . In the case of the form factors this is not yet possible. 

The procedure employed in this paper is to draw on the 
earlier phenomenological experience with chiral extrapo- 
lations of the magnetic moments and charge radii, taking 
simple phenomenological functional forms which build in 
the correct non-analytic behaviour and any other asymp- 
totic constraints that are known |29,. For the present the 
accuracy of the lattice data is not such that it could dis- 
criminate between different functional forms for the Q^- 
depcndcnce of the nucleon form factors. In particular, it 
cannot yet address the JLab issue of whether Ge/Gm 
decreases as increases. We therefore parametrize sep- 
arately the isoscalar and isovector lattice data, at a given 
value of TO„, or m^, as a dipole and then extrapolate the 
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dipole mass as a function of m^r , building in the appropri- 
ate chiral constraints. Clearly this simple approach can be 
systematically improved as the lattice simulations become 
capable of making a better discrimination between possi- 
ble functional forms. For the present we shall see that our 
relatively simple approach already produces quite impres- 
sive results. 



2 Extrapolations 

It is well known ^ that the Q^-dependence of the nucleon 
Sachs electromagnetic form factors (with the exception of 
the neutron electric form factor G^) are described in first 
approximation by a dipole form 



G(0) 



(1) 



Here A is the dipole mass and the charge form factor of 
the proton satisfies G^(0) = 1, while G^(0) = fip and 
G^(0) = Hn are the proton and neutron magnetic mo- 
ments. As explained in the Introduction, we will use this 
phenomenological fact to construct a simple but effective 
extrapolation formula. 

To isolate the chiral behaviour of the form factors we 
rearrange them into isovector and isoscalar combinations: 



G" = GP - G" (2) 



and 



G' = GP + G", (3) 

respectively. These also display a dipole- like Q^-dependence 
but with different dipole masses and magnetic moments. 

To extrapolate lattice QCD results for the electromag- 
netic form factors from the large pion masses at which they 
are calculated to the physical regime, we extract dipole 
masses and magnetic moments from the lattice data and 
then extrapolate these as a function of m^. 

Following Refs. |11^HT^ . one can use a Fade approxi- 
mant which builds in both the correct chiral non-analytic 
behaviour as to^ and the correct asymptotic be- 
haviour as rriq — > oo to extrapolate the neutron and proton 
magnetic moments 



Mo 



1 



(4) 



The chiral coefficients for the isovector and isoscalar mo- 
ments are Xv = —8.82 and Xs = respectively and /io and 
c are fitting parameters, to be determined by the lattice 
data. 

In order to build a suitable extrapolating function for 
the Q^-dependence of the form factors, we use the connec- 
tion between the mass parameter in a dipole form factor 
and the corresponding mean-square radius. For the isovec- 
tor magnetic form factor the mean square radius is 
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Fig. 1. Fit to values of the isovector electric form factor 
dipole mass extracted from lattice data with lattice spacing 
a = 0.051 fm. The physical value predicted by the fit is also 
indicated. 
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Fig. 2. Fit to values of the isovector magnetic form factor 
dipole mass extracted from lattice data with lattice spacing 
a = 0.051 fm. The physical value predicted by the fit is also 
indicated. 



Comparing this to the dipole of Eq. we find an ex- 
pression relating the dipole mass to {r'^) 



V \2 

m) 



12 



(6) 



M 



The chiral behaviour of the magnetic mean squared radius 
is known from chiral perturbation theory j3U| 



Xi , , ,m^. 
\-X2 ln( ). 



The constants xi and X2 are given by 

Xi 



9a^n 



87r/2K„ ' 



(7) 
(8) 
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where = 1.27 is the axial coupUng constant and = 
93 MeV is the pion decay constant, niN = 940 MeV is 
the nucleon mass and k.^ ~ 4.2 is the isovector anomalous 
magnetic moment of the nucleon (in the chiral limit). 
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Fig. 3. Linear fit to values of the isoscalar electric form factor 
dipole mass extracted from lattice data with lattice spacing 
a = 0.051 fm. The physical value predicted by the fit is also 
indicated. 
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Fig. 4. Linear fit to values of the isoscalar magnetic form factor 
dipole mass extracted from lattice data with lattice spacing 
a = 0.051 fm. The physical value predicted by the fit is also 
indicated. 



Earlier experience with nucleon properties as a func- 
tion of quark mass has suggested that, as a consequence 
of the finite size of the source of the pion field, chiral loops 
are strongly suppressed for > 0.4 GeV [HJ- Thus in 
order to build an effective chiral extrapolation formula 
one needs to modify the non-analytic terms given above 
so that they are suppressed above this mass. For guidance 
as to an appropriate parametrization of this suppression 
we have evaluated the pion loops that give rise to the non- 



analytic chiral terms using the Cloudy Bag Model (CBM) 
tJ7:'3J'^~4'. The results led us to replace Eq. 10 with the 
expression 

i^Xi - — -arctan(Ai/m.) + ^ln(^^|-^), (10) 

which ensures the correct chiral behaviour at low-m^r but 
suppresses it at values larger than /j,. 

6.0 I \ 1 \ 1 \ \ \ 1 
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Fig. 5. Fit to values of the isovector magnetic moment ex- 
tracted from lattice data with lattice spacing a — 0.051 fm. 
The physical value predicted by the fit is also indicated. 
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Fig. 6. Fit to values of the isoscalar magnetic moment ex- 
tracted from lattice data with lattice spacing a = 0.051 fm. 
The physical value predicted by the fit is also indicated. 



Substituting this functional form into Eq. 10) and in- 
troducing a linear dependence on to give the expected 
behaviour at large quark mass, we get the following ex- 
pression for the dipole mass associated with the isovector 
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magnetic form factor 



[All? 



12(l + ylim2) 



^0 



XI 2 



— arctan(/i/mTr) 



2 



(11) 

Here Aq and Ai are unknown parameters, adjusted to fit 
the lattice data. 

We can perform a similar analysis to find an expres- 
sion for the isovector electric form factor dipole mass A^^ , 
combining the non-analytic chiral behaviour predicted by 
xPT at low with the slowly varying behaviour ob- 
served in lattice QCD at large masses. Using the equiva- 
lent of Eq. © for the electric form factor, we can relate 
the dipole mass to the mean charge radius by 



12 



(12) 



The chiral behaviour of the isovector mean charge radius 
is like that of the magnetic radius but with no 1/ ttItt term 



(13) 



Combining this chiral behavior with a linear dependence 
on produces the following predicted form for the isovec- 
tor electric form factor dipole mass, 



12{l + Biml) 



Bn 



X2 



) 



(14) 



where Bq and Bi are general fitting parameters. 

Having isolated the chiral non-analytic behaviour to 
the isovector form factors, we expect that the isoscalar 
dipole masses should be roughly linear in - as observed 
in CBM results. 
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Fig. 7. The magnetic form factor for the proton extrapolated 
from lattice data with lattice spacings a = 0.093 fm (single 
line), a = 0.068 fm (dash-dot line) and a — 0.051 fm (dotted 
line) and compared to experimental results. 
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Fig. 8. The magnetic form factor for the neutron extrapolated 
from lattice data with lattice spacings a = 0.093 fm (single 
line), a — 0.068 fm (dash-dot line) and a — 0.051 fm (dotted 
line) and compared to experimental results. 



The forms chosen to represent the dipole masses as a 
function of are able to reproduce the predictions of 
the CBM El very weU. 



3 Results 

The QCDSF Collaboration T has recently reported QQCD 
results for the isovector and isoscalar nucleon electric and 
magnetic form factors calculated at different values 
and pion masses and for three different lattice spacings 
(/3 = 6.0, 6.2 and 6.4). Using this data, with the scale set 
using the Sommer method (ro = 0.5fm [SS|)j we have plot- 
ted the form factors G^j, G%j, and G%, calculated on 
the lattice as a function of Q^, at each pion mass and lat- 
tice spacing, and fit each graph with a dipole form, finding 
a best-fit dipole mass and magnetic moment. Plots of these 
best-fit dipole masses and magnetic moments against 
for the lattice spacing a — 0.051 fm are shown in Figs.^ 
HElandll 

The electric and magnetic isovector dipole masses are 
fitted with functions of the form Eq. (|ll|l and Eq. H14|> . re- 
spectively. The scale fi = 0.41 GeV was chosen to give the 
best simultaneous fit to the lattice data for both the elec- 
tric and magnetic isovector form factors at all three lattice 
spacings. We observe that the curves fit the data well and 
predict physical results significantly different from those 
obtained by a naive linear extrapolation. In both cases, the 
form of the fitting functions greatly restricts the physical 
value, resulting in a small error. Finally, we note that the 
electric and magnetic isoscalar dipole mass plots are fitted 
with linear functions in m^, because of the absence of any 
leading non-analytic behaviour in this case. 

The magnetic moment plots are extrapolated using 
the Fade approximant, Eq. Q. The resulting errors in 
the physical magnetic moments dominate over errors in 
the dipole mass and are the main source of error in our 
results for the physical magnetic form factors. Similar 
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Fig. 9. The electric form factor for the proton extrapolated 
from lattice data with lattice spacings a = 0.093 fm (single 
line), a = 0.068 fm (dash-dot line) and a — 0.051 fm (dotted 
line) and compared to experimental results. 
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Fig. 10. The electric form factor for the neutron extrapolated 
from lattice data with lattice spacings a = 0.093 fm (single 
line), a = 0.068 fm (dash-dot line) and a — 0.051 fm (dotted 
line) compared to experimental results. 



fits were produced for lattice results at lattice spacings 
a = 0.068 fm and a = 0.093 fm. 

Using the extrapolated physical values for the dipole 
masses and magnetic moments, we reconstruct the electric 
and magnetic form factors for the proton and neutron as 
a sum of two dipoles in Q^. Figs. [3 |S1 151 and ^| show 
the electric and magnetic form factors for the proton and 
neutron extrapolated from that lattice data at each of the 
three lattice spacings (a = 0.093 fm, a = 0.068 fm and 
a = 0.051 fm) compared to the experimental data from 
Refs. ISni. 

Except for the case of the neutron electric form factor 
(which is very sensitive because it is the result of a cancel- 
lation between and G|; which yields a relatively small 
result) our extrapolated curves match the experimental 
points very well, particularly for and G^j. Although 



our curve for is almost twice as big as the experimen- 
tal values, it does describe the correct shape but peaks at 
slightly too high a value of Q^. 



4 Conclusion 

We have proposed a relatively simple approach to the ex- 
trapolation of lattice QCD data for the nucleon electro- 
magnetic form factors. The data from QCDSF has been 
parametrized by a simple dipole form, with the dipole 
mass parameter taken to be a function of ttItt which incor- 
porates the leading non-analytic behaviour of chiral per- 
turbation theory. For the isoscalar case, where there is no 
leading non-analytic behaviour, the extrapolation is a sim- 
ple linear function of m^, while for the isovector case we 
have used a functional form suggested by studies based on 
the cloudy bag model, which suppress the chiral behaviour 
for pion masses larger than 400 MeV. 

The level of agreement between the empirical de- 
pendence of the proton electric form factor and the proton 
and neutron magnetic form factors is impressive. This is 
quite remarkable when one considers that the data is based 
on quenched approximation and is extrapolated from a 
light quark mass a factor of twenty above the physical 
mass. This is only possible because of the remarkable fact 
that pion loops are suppressed by finite size effects once 
the pion Compton wavelength is smaller than the size of 
the source - empirically, 111^^ > 0.4 GeV. The extrapola- 
tion of the data has, of course, been performed using chiral 
coefficients appropriate to full QCD. 

One of the main open questions in the analysis con- 
cerns the neutron electric form factor, which will always 
be a challenge given that it vanishes at = - as a re- 
sult of the cancellation between the isoscalar and isovector 
contributions. In addition, we observe that there appears 
to be some residual (3 dependence in the value of the nu- 
cleon magnetic moments and this is the major ambiguity 
in the current reconstruction of the proton and neutron 
magnetic form factors. It is encouraging that the results 
for P = 6.4, corresponding to the smallest lattice spac- 
ing, is in the best agreement with the experimental data. 
Future lattice simulations will undoubtedly clarify this is- 
sue and also give us data at lower values of the quark 
mass. It would also be valuable to have full QCD sim- 
ulations (rather than QQCD) to ensure that there is no 
unanticipated systematic difference. Most importantly, as 
the accuracy of the data and the range of covered in- 
creases, we can extend the present analysis by employing 
more complicated fitting functions which should allow us 
to test the behaviour of properties such as Ge/ Gm against 
experiment. 
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